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Abstract 

In this paper, we study the joint resource allocation algorithm design for downlink and uplink 
multicarrier transmission assisted by a shared user equipment (UE)-side distributed antenna system 
(SUDAS). The proposed SUDAS simultaneously utilizes licensed frequency bands and unlicensed 
frequency bands, (e.g. millimeter wave bands), to enable a spatial multiplexing gain for single-antenna 
UEs to improve energy efficiency and system throughput of 5-th generation (5G) outdoor-to-indoor 
communication. The design of the UE selection, the time allocation to uplink and downlink, and the 
transceiver processing matrix is formulated as a non-convex optimization problem for the maximization 
of the end-to-end system energy efficiency (bits/Joule). The proposed problem formulation takes into 
account minimum data rate requirements for delay sensitive UEs and the circuit power consumption 
of all transceivers. In order to design a tractable resource allocation algorithm, we first show that the 
optimal transmitter precoding and receiver post-processing matrices jointly diagonalize the end-to-end 
communication channel for both downlink and uplink communication via SUDAS. Subsequently, the 
matrix optimization problem is converted to an equivalent scalar optimization problem for multiple 
parallel channels, which is solved by an asymptotically globally optimal iterative algorithm. Besides, 
we propose a suboptimal algorithm which finds a locally optimal solution of the non-convex optimization 
problem. Simulation results illustrate that the proposed resource allocation algorithms for SUDAS 
achieve a significant performance gain in terms of system energy efficiency and spectral efficiency 
compared to conventional baseline systems by offering multiple parallel data streams for single-antenna 
UEs. In fact, the proposed SUDAS is able to bridge the gap between the current technology and the 
high data rate and energy efficiency requirements of 5G outdoor-to-indoor communication systems. 
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I. Introduction 

High data rate, high energy effieiency, and ubiquity are basic requirements for 5-th generation 
(5G) wireless communication systems. A relevant technique for improving the system throughput 
for given quality-of-service (QoS) requirements is multiple-input multiple-output (MIMO) IfTI- 
a, as it provides extra degrees of freedom in the spatial domain which facilitates a trade-off 
between multiplexing gain and diversity gain. In particular, massive MIMO, which equips the 
transmitter with a very large number of antennas to serve a comparatively small number of user 
equipments (UEs), has received considerable interest recently [|2l, [|3l. The high flexibility in 
resource allocation makes massive MIMO a strong candidate for 5G communication systems. 
However, state-of-the-art UEs are typically equipped with a small number of receive antennas 
which limits the spatial multiplexing gain offered by MIMO to individual UEs. On the other 
hand, the combination of millimeter wave (mmW) and small cells, e.g. femtocells, has been 
proposed as a core network architecture for 5G indoor communication systems flU, [|51 since 
most mobile data traffic is consumed indoors [l6]|. The huge free, unlicensed frequency spectrum 
in the mmW frequency bands appears to be suitable and attractive for providing high speed 
communication services over short distances in the order of meters. However, the backhauling 
of the data from the service providers to the small cell base stations (BSs) is a fundamental system 
bottleneck. In general, the last mile connection from a backbone network to the UEs at homes 
can only support high data rates if optical fibers are deployed, which is known as fiber-to-the- 
home (ETTH). Yet, the cost in deploying ETTH for all indoor users is prohibitive. Eor instance, 
the cost in equipping every building with ETTH in Germany is estimated to be around 67 billion 
Euros ATI which makes high speed small cells not an appealing universal solution for 5G indoor 
communication systems in terms of implementation cost. Another attractive system architecture 
for 5G is to combine massive MIMO with mmW communications AH, All by using outdoor 
mmW BSs. However, the high penetration loss of building walls limits the suitability of mmW 
for outdoor-to-indoor communication scenarios. Thus, additional effective system architecture 
for outdoor-to-indoor communication is needed. 

Distributed antenna systems (DAS) are an existing system architecture on the network side 
and a special form of MIMO. DAS are able to cover the dead spots in wireless networks, extend 
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service coverage, improve spectral efficiency, and mitigate interference IfTOll . IfTTlI . It is expected 
that DAS will play an important role in 5G communication systems ifT^ . Specifically, DAS can 
realize the potential performance gains of MIMO systems by sharing antennas across different 
terminals of a communication system to form a virtual MIMO system IfTSlI . Lately, there has 
been a growing interest in combining orthogonal frequency division multiple access (OFDMA) 
and DAS to pave the way for the transition of existing communication systems to 5G ifTdl - lfT^ . 
In lfT4l . the authors studied suboptimal resource allocation algorithms for multiuser MIMO- 
OFDMA systems. In ifTSll . a utility-based low complexity scheduling scheme was proposed 
for multiuser MIMO-OFDMA systems to strike a balance between system throughput and 
computational complexity. The optimal subcarrier allocation, power allocation, and bit loading 
for OFDMA-DAS was investigated in [fT^ . However, similar to massive MIMO, DAS cannot 
significantly improve the data rate of individual UE when the UEs are single-antenna devices. 
Besides, the results in lfT4l - in^ . which are valid for either downlink or uplink communication, 
may no longer be applicable when joint optimization of downlink and uplink resource usage 
is considered. Eurthermore, the total system throughput in [[T4ll - [[T^ is not only limited by 
the number of antennas equipped at the individual UEs, but is also constrained by the system 
bandwidth which is a very scarce resource in licensed frequency bands. In fact, licensed spectrum 
is usually located at sub 6 GHz frequencies which are suitable for long distance communication. 
On the contrary, the unlicensed frequency spectrum around 60 GHz offers a large bandwidth 
of 7 GHz for wireless communications but is more suitable for short distance communication. 
The simultaneous utilization of both licensed and unlicensed frequency bands for high rate 
communication introduces a paradigm shift in system and resource allocation algorithm design 
due to the related new challenges and opportunities. Yet, the potential system throughput gains of 
such hybrid systems have not been thoroughly investigated in the literature. Thus, in this work, we 
study the resource allocation design for hybrid communication systems simultaneously utilizing 
licensed and unlicensed frequency bands to improve the system performance. 

An important requirement for 5G systems is energy efficiency. Over the past decades, the 
development of wireless communication networks worldwide has triggered an exponential growth 
in the number of wireless communication devices for real time video teleconferencing, online 
high definition video streaming, environmental monitoring, and safety management. It is expected 
that by 2020, the number of interconnected devices on the planet may reach up to 50 billion 
uni. The related tremendous increase in the number of wireless communication transmitters 
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and receivers has not only led to a huge demand for licensed bandwidth but also for energy. 
In particular, the escalating energy consumption of electronic circuitries for communication and 
radio frequency (RF) transmission increases the operation cost of service providers and raises 
serious environmental concerns due to the produced green house gases. As a result, energy 
efficiency has become as important as spectral efficiency for evaluation of the performance of 
the resource utilization in communication networks. A tremendous number of green resource 
allocation algorithm designs have been proposed in the literature for maximization of the energy 
efficiency of wireless communication systems Q, ifT^ - lEOll . In Q, joint power allocation and 
subcarrier allocation was considered for energy-efficient massive MIMO systems. In [fT^ . the 
energy efficiency of a three-node multiuser MIMO system was studied for the two-hop compress- 
and-forward relaying protocol. The trade-off between energy efficiency and spectral efficiency 
in DAS for fair resource allocation in flat fading channels was studied in ifT^ . Power allocation 
for energy-efficient DAS was investigated in ll20ll for frequency-selective channels. However, in 
0, lfT^ - ll20ll . it was assumed that the transmit antennas were deployed by service providers and 
are connected to a central unit by high cost optical fibers or cables for facilitating simultaneous 
transmission which may not be feasible in practice. To avoid this problem, unlicensed and 
licensed frequency bands may be used simultaneously to create a wireless data pipeline for 
DAS to provide high rate communication services. Nevertheless, the resource allocation algorithm 
design for such a system architecture has not been investigated in the literature, yet. 

In this paper, we propose a shared UE-side distributed antenna system (SUDAS) to assist 
the outdoor-to-indoor communication in 5G wireless communication systems. In particular, 
SUDAS simultaneously utilizes licensed and unlicensed frequency bands to facilitate a spatial 
multiplexing gain for single-antenna transceivers. We formulate the resource allocation algorithm 
design for SUDAS assisted OFDMA downlink/uplink transmission systems as a non-convex 
optimization problem. By exploiting the structure of the optimal precoding and post-processing 
matrices adopted at the BS and the SUDAS, the considered matrix optimization problem is 
transformed into an equivalent optimization problem with scalar optimization variables. Capital¬ 
izing on this transformation, we develop an iterative algorithm which achieves the asymptotically 
globally optimal performance of the proposed SUDAS for high signal-to-noise ratios (SNRs) and 
large numbers of subcarriers. Also, the asymptotically optimal algorithm serves as a building 
block for the design of a suboptimal resource allocation algorithm which achieves a locally 
optimal solution for the considered problem for arbitrary SNRs. 
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Fig. 1. The upper half of the figure illustrates the downlink and uplink communication between a base station (BS) and K = 3 
user equipments (UEs) assisted by M = 3 SUDACs. The proposed system utilizes a licensed frequency band and an unlicensed 
frequency band such as the mmW band (e.g. ~ 60 GHz). The lower half of the figure depicts the time division duplex (TDD) 
approach adopted for downlink and uplink communication within a coherence time slot. 

II. SUDAS Assisted OFDMA Network Model 
A. SUDAS System Model 

We consider a SUDAS assisted OFDMA downlink (DL) and uplink (UL) transmission network 
which consists of one N antenna BS, a SUDAS, and K single-antenna UEs, cf. Figure [B The 
BS is half-duplex and equipped with N antennas for transmitting and receiving signals in a 
licensed frequency band. The UEs are single-antenna devices receiving and transmitting signals 
in the unlicensed frequency band. Also, we focus on a wideband multicarrier communication 
system with np orthogonal subcarriers. A SUDAS comprises M shared user equipment (UE)-side 
distributed antenna components (SUDACs). A SUDAC is a small and cheap device deployed 
inside a building^ which simultaneously utilizes both a licensed and an unlicensed frequency band 
for increasing the DL and UL end-to-end communication data rate. A basic SUDAC is equipped 
with one antenna for use in a licensed band and one antenna for use in an unlicensed band. 
We note that the considered single-antenna model for SUDAC can be extended to the case of 
antenna arrays at the expense of a higher complexity and a more involved notation. Eurthermore, 
a SUDAC is equipped with a mixer to perform frequency up-conversion/down-conversion. Eor 
example, for DL communication, the SUDAC receives the signal from the BS in a licensed 
frequency band, e.g. at 800 MHz, processes the received signal, and forwards the signal to the 


*In practice, a SUDAC could be integrated into electrical devices such as electrical wall outlets, switches, and light outlets. 
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Fig. 2. Illustration of signal forwarding from(/to) a licensed band to(/from) different unlicensed frequency sub-bands in the 
SUDAS. 

UEs in an unlicensed frequeney band, e.g. the mmW bands. We note that sinee the BS-SUDAC 
link operates in a sub-6 GHz lieensed frequeney band, it is expeeted that the assoeiated path loss 
due to bloekage by building walls is mueh smaller eompared to the ease where mmW bands were 
direetly used for outdoor-to-indoor eommunieation. Henee, the BS-to-SUDAS ehannel serves as 
wireless data pipeline for the SUDAS-to-UE eommunieation ehannel. Also, sinee signal reeeption 
and transmission at eaeh SUDAC are separated in frequeney, ef. Eigure[2]and ETI . simultaneous 
signal reeeption and transmission ean be performed in the proposed SUDAS whieh is not possible 
for traditional relaying systemj^ due to the limited speetrum availability in the lieensed bands. 
The UL transmission via SUDAS ean be performed in a similar manner as the DL transmission 
and the detailed operation will be diseussed in the next seetion. In praetiee, a huge bandwidth 
is available in the unlieensed bands. Eor instanee, there is nearly 7 GHz unlieensed frequeney 
speetrum available for information transmission in the 57 — 64 GHz band (mmW bands). In 
this paper, we study the potential system performanee gains for outdoor-to-indoor transmission 
aehieved by the proposed SUDAS arehiteeture. In partieular, we foeus on the ease where the 
SUDACs are installed in eleetrieal wall outlets indoor and ean eooperate with eaeh other by 
sharing ehannel state information, power, and reeeived signals, e.g. via power line eommunieation 
links. In other words, for the proposed resouree alloeation algorithm, joint proeessing aeross the 
SUDACs is assumed to be possible sueh that the SUDACs ean fully exploit the degrees of 
freedom offered by their antennas. The joint proeessing arehiteeture of the SUDAS in this paper 


^Since the BS-to-SUDAS and SUDAS-to-UE links operate in two different frequency bands, the proposed SUDAS should 
not be considered a traditional relaying system (22). 
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reveals the maximum potential performance gain of the proposed SUDAS. 

Furthermore, we adopt time division duplexing (TDD) to facilitate UL and DL communication 
for half-duplex UEs and BS. To simplify the following presentation, we assume a normalized unit 
length time slot whose duration is the coherence time of the channel, i.e., the communication 
channel is time-invariant within a time slot. Each time slot is divided into two intervals of 
duration, a and which are allocated for the DL and UL communication, respectively. 


B. SUDAS DL Channel Model 


In the DL transmission period a, the BS performs spatial multiplexing in the licensed band. 
The data symbol vector G on subcarrier z G {1, ..., np} for UE /c G {1, ..., K} 

is precoded at the BS as 


X 


[i,k] 

DL 


■^DL ^DL ) 


( 1 ) 


where G is the precoding matrix adopted by the BS on subcarrier i and C^^^® 

denotes the set of all N x matrices with complex entries. The signals received on subcarrier 
i at the M SUDACs for UE k are given by 


[i,k\ 

Ts-dl 


xq-m I 

-^B^S-^DL ^ 


( 2 ) 


where = [2/s-dli> ■ ■ ■ >?/s-dlm]^’ ?/s-dl„, denotes the received signal at SUDAC m G 

M}, and (•)^ is the transpose operation. Hg^g is the M x N MIMO channel matrix 
between the BS and the M SUDACs on subcarrier i and captures the joint effects of path 
loss, shadowing, and multi-path fading, is the additive white Gaussian noise (AWGN) vector 
impairing the M SUDACs in the licensed band on subcarrier i and has a circularly symmetric 
complex Gaussian (CSCG) distribution CAA(0, S) on subcarrier i, where 0 is the mean vector 
and S is the M x M covariance matrix which is a diagonal matrix with each main diagonal 
element given by Nq. 

In the unlicensed band, each SUDAC performs orthogonal/re^wency repetition. In particular, 
the M SUDACs multiply the received signal vector on subcarrier i, ygi^^L, by G and 

forward the processed signal vector to UE k on subcarrier z in M different independent frequency 
sub-bands in the unlicensed spectrurry, cf. Eigure |2l In other words, different SUDACs forward 


'^For a signal bandwidth of 20 MHz, there can be 350 orthogonal sub-bands available within 7 GHz of bandwidth in the 
60 GHz mmW band llH. For simplicity, we assume that each of the M SUDACs uses one fixed sub-band for DL and UL 
communication. 
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that received signals in different sub-bands and thereby avoid multiple access interference in the 
unlicensed spectrum. 

The signal received at UE k on subcarrier i from the SUDACs in the M frequency bands, 
Ys^ue ^ , can be expressed as 


[i,k] 

Ys-^ue 




pfijfe] I TTpifc] 

^DL 


S^UE-^ DL 




desired signal 


amplified noise 


(3) 


\i k] 

The m-th element of vector YsVue represents the received DL signal at UE k in the m- th 
unlicensed frequency sub-band. Since the SUDACs forward the DL received signals in different 
orthogonal frequency bands, is a diagonal matrix with the diagonal elements representing 

the channel gain between the SUDACs and UE k on subcarrier i in unlicensed sub-band m. 

G is the AWGN vector at UE k on subcarrier i with distribution CAA(0, S^), where 

Sfc is an M X M diagonal matrix and each main diagonal element is equal to A^uEfe- 
We assume that M > and UE k employs a linear receiver for estimating the DL data vector 
symbol received in the M different sub-bands in the unlicensed band. Hence, the estimated data 
vector symbol, G on subcarrier i at UE k is given by 



(wK')" 


[i,k] 

Ys^-uei 


(4) 


where G is a post-processing matrix used for subcarrier i at UE k, and (■)^ denotes 

the Hermitian transpose. Without loss of generality, we assume that (dpL^)^} = Ia^s 

where Ijvg is an Ng x Ng identity matrix and £{■} denotes statistical expectation. As a result, 
the minimum mean square error (MMSE) matrix for data transmission on subcarrier i for UE k 
via the proposed SUDAS and the optimal MMSE post-processing matrix are given by 

_ fffpi[i,k] _ ^[i,k]\(AliM _ ^li,k]\Hx _ 


-'DL ~ ^U'-‘dL '-‘dL If-^DL "-^DL 


and W&JI = + 0K')-‘r&?, 


J - ^Ns + DL > 1'='dL ) ^ DL 


-1 


(5) 

( 6 ) 




\i k] 

denotes the matrix inverse, ECy is the effective end-to-end channel 


respectively, where 

matrix from the BS to UE k via the SUDAS on subcarrier i, and is the corresponding 
equivalent noise covariance matrix. These matrices are given by 

H 


[*i^] _ TT[*,fc] 17 [*ifc]Tj[*] 


i,k] 


DL 


TTPAJ Tp[‘AjTT[»J p 

■^s->-ue-^dl -^b^-s-^dl 


nnH — I pl 

ana 


■[i,k] Tn[i,A:] 


■^S^-UE-^ DL 


+ Im. (7) 
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Remark 1: The SUDAS concept is fundamentally different from traditional relaying systems 
which aim at extending service coverage lf2^ . For DL communication, the SUDAS converts 
the spatial multiplexing performed at the BS in the licensed band into frequency multiplexing 
in the unlicensed band to allow single-antenna UEs to decode multiple spatial data streams. 


C. SUDAS UL Channel Model 


In the UL transmission period (3, UE k performs frequency multiplexing in the unlicensed 
band. The data symbol vector G on subcarrier i G {1, ..., np} from UE k is 

precoded as 



^UL 


^UL ; 


( 8 ) 


where G is the UL precoding matrix adopted by UE k on subcarrier i over the M 

different frequency sub-bands in the unlicensed spectrum. The signals received on subcarrier i 
at the M SUDACs for UE k are given by 

yS-UL “ ■'^UE-)-S^UL + 

where = [?/s3jli • • • 2/s-ulm]^’ ys-vh^ denotes the received signal at SUDAC m in 

unlicensed frequency sub-band m G M}, and is the AWGN impairing the M 

SUDACs on subcarrier i in the unlicensed frequency band, has distribution CAf{0, Sul), 
where Sul is an M x M diagonal matrix and each main diagonal element is equal to Aul- 
H&eLs is a diagonal matrix with the main diagonal elements representing the channel gains 
between UE k and the M SUDACs on subcarrier i in unlicensed sub-band m. In fact, the UEs- 
to-SUDAS channels serve as a short distance wireless data pipeline for the SUDAS-to-BS UL 
communication. 

Each SUDAC forwards the signals received in the unlicensed band in the licensed band to 
assist the UL communication. In particular, the M SUDACs multiply the received signal vector 
on subcarrier i by G and forward the processed signal vector to the BS on subcarrier 

i in the licensed spectrum, cf. Eigure |2l As a result, the signal received at the BS from UE k 
on subcarrier i via the SUDAS, ^ can be expressed as 


[i,k] 

y^B ^ 


UL ■‘^UE^-S-'^UL "-^UL 

'-V-' 

desired signal 


i,k\ ^ 


F 

■^S^>B-^ UL 


+nl5’"l 


( 10 ) 


amplified noise 


U1 U ^1 

Matrix is the UL channel between the M SUDACs and the BS on subcarrier i, and rig ‘ 

is the AWGN vector in subcarrier i at the BS with distribution CAA(0, XIb), where Sb is an 
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M X M diagonal matrix and each main diagonal element is equal to A^b- At the BS, we assume 
that N > Ns and the BS employs a linear receiver for estimating the data vector symbol received 
from the SUDAS in the licensed band. The estimated data vector symbol, cIul' G on 

subcarrier i at the BS from UE k is given by 



(W&a- 




(11) 


where G is a post-processing matrix used for subcarrier i at UE k. Without loss 

of generality, we assume that £^{duL^(duL')^} = Iats- As a result, the MMSE matrix for data 
transmission on subcarrier i from UE k to the BS via the SUDAS and the optimal MMSE 
post-processing matrix are given by 


eK' = £{(<i 


i,k] 

UL 


- ^UL KaUL - “UL ) J - -'-Ars + (i UL . 


and = 


A UL U UL 1 '-^UL 




UL 


i,k] 

UL 


-1 


Lfc]\-l-p[Lfc] 


UL ) 


( 12 ) 

(13) 


H k] 

respectively, where is the effective end-to-end channel matrix from UE k to the BS via 

\i k] 

the SUDAS on subcarrier z, and ©ul is the corresponding equivalent noise covariance matrix. 
These matrices are given by 


^ UL 


ttW td 

UL ■'^UE-j-S-'^^UL 


i,k] 


and © 


i,k] 

UL 


TT[*,fc] 

■•^UE-j-S-*^ UL 


TtIL^] Tp[*ifc] 
■■^UE^-S-*^ UL 


H 


+ Im-(14) 


Remark 2: Since TDD is adopted and DE and UE transmission occur consecutively within 
the same coherence time, for resource allocation algorithm design, it is reasonable to assume 
that channel reciprocity holds, i.e., = (Hg^^g)^ and = (Hg^g)^. 

III. Problem Eormulation 

In this section, we first introduce the adopted system performance measure. Then, the design 
of resource allocation and scheduling is formulated as an optimization problem. 

A. System Throughput 

The end-to-end DE and UE achievable data rate on subcarrier i between the BS and UE k 
via the SUDAS are given by |[25]l 

^&L = - log 2 (det[E£’L']) and = - loga (det[ES;f']) , (15) 
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respectively, where det(-) is the determinant operation. The DL and UL data rate (bits/s) for UE 
k can be expressed as 



*DL 


r)[i,k] 

-'^UL 


and pPl = Y 1 


np 


'’UL -'^DL ) 


(16) 


2=1 2=1 

respectively, where e {0,a} and e {0,/3} are the discrete subcarrier allocation 
indicators, respectively. In particular, a DL and an UL subcarrier can only be utilized for a 
and (3 portions of the coherence time, respectively, or not be used at all. 

The system throughput is given by 


i'P, *5) = dSL + pPl [bits/s], 

k=l k=l 


(17) 


where V = {P^’l , Pul', Ful'} and S = {sul', Sul', a, /?} are the precoding and subcarrier 
allocation policies, respectively. 

On the other hand, the power consumption of the considered SUDAS assisted communication 
system consists of seven terms which can be divided into three groups and expressed as 


Wtp(P,5) - Pcb + ^-PAntB + dTPcsuDAC + 


System circuit power consumption 

K n-p K Up 


EE4LeBTr(pM(pM)« 

k=l i=l 


k=l i=l 


Total DL transmit power consumption 

K riY K np 

E c/. E 4 'l tv (pKipliL')") + E E iv ( ga 


k=l i=l 


k=l i=l 


Total UL transmit power consumption 


(18a) 

(18b) 

[Joule/s] (18c) 


where 


DL 


hLsPdl(Pdl)''(hLs)" 


l-M 


A DL ) ) 


(19) 


G 


'i,k] 

UL 


= F 


[i,k] /Trfijfe] 


UL 


■^UE^.S-*^^UL 1-*^^UL 1 


H, 


H 


[i,k] 
VE^SJ 


\H 


•-M 


A UL ) ) 


( 20 ) 


and Tr(-) is the trace operator. The three positive constant terms in (I18al) . i.e., Pcb^ Pcsvoac^ 
and PcuE’ represent the power dissipation of the circuits lf2^ for the basic operation of the BS, 
the SUDAC, and the UE, respectively, and PAntB denotes the circuit power consumption per BS 
antenna. Equations (I18bl) and (I18cl) denote the total DL transmit power consumption and the 
total UL power consumption, respectively. Specifically, Tr(GuL') and Tr(GuL') are the DL and 
UL transmit powers of the SUDAS needed for facilitating the DL and UL communication of UE 
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k in subcarrier i, respectively. Similarly, Tr ( Pdl^(Pdl^)^ ) and Tr ( Pul'(Pul')^ ) represent 


the DL transmit power from the BS to the SUDAS for UE k and the UL transmit power from 
UE k to the SUDAS in subcarrier i, respectively. To capture the power inefficiency of power 
amplifiers, we introduce linear multiplicative constants £b,£S) and £k for the power radiated by 
the BS, the SUDAS, and UE k in (fTSl) . respectively. Eor instance, if £b = 4, then for 1 Watt of 
power radiated in the RE, the BS consumes 4 Watt of power which leads to a power amplifier 
efficiency of 259^^. 


The energy efficiency of the considered system is defined as the total number of bits exchanged 
between the BS and the K UEs via the SUDAS per Joule consumed energy: 


UesiV,S) = [bits/Joule]. 


( 21 ) 


B. Problem Formulation 

The optimal precoding matrices, V* = {P^l^*, and the optimal sub¬ 
carrier allocation policy, S* = a*,/?*}, can be obtained by solving the following 

optimization problem: 

maximize Uesi'P.S) 

'P,S 


K np 


.t. Cl: 4 ’l Tr I 


k=l i=l 
K np 


DL 


< Ft, 


EE 


hk] ^ qM] I < 


’DL 


k=l i=l 

np 


max 5 


C3: 4^1 IV (pSIpK')") < P„„,Vfc e {1,..., K], 


2=1 
K np 


«: E E TV (gK') < P™ , 


k=l i=l 
Sk\ 


C5: C6:4[>P™.,V)!6® 


np 


C7: yy4v‘><a.Vi, 


np 


2=1 


C9: 4^'e{0,a},V^,A:, 

Cll: « + /?<!, 


C8:5^4’l'’</3, V^, 

2=1 

CIO: 4J’e{0,/3},Vz,A:, 
C12: a,/? > 0. 


( 22 ) 


this paper, we assume that Class A power amplifiers with linear characteristic are implemented at the transceivers. The 
maximum power efficiency of Class A amplifiers is limited to 25%. 
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Constants Pt and in Cl and C2 are the maximum transmit power allowanees for the 

BS and the SUDAS (M SUDACs) for DL transmission, respeetively, where P^ax is the average 
transmit power budget for a SUDAC. Similarly, eonstraints C3 and C4 limit the transmit power 
for UE k and the SUDAS (M SUDACs) for UL transmission, respeetively, where Pmax^ and 
are the maximum transmit power budgets of UE k and the SUDAS, respeetively. We note that 
in praetiee the maximum transmit power allowanees for the SUDAS-to-UE, Pmax> and SUDAS- 
to-BS, Pmax’ ^aiay be different due to different regulations in lieensed and unlieensed bands. Sets 
Pdl and Pul in eonstraints C5 and C6 denote the set of delay sensitive UEs for DE and UE 
eommunieation, respeetively. In partieular, the system has to guarantee a minimum required DE 
data rate and UE data rate if UE k requests delay sensitive serviees in the DE 

and UE, respeetively. Constraints C7 - CIO are imposed to guarantee that eaeh subearrier ean 
serve at most one UE for DE and UE eommunieation for fraetions of a and (3 of the available 
time. Constraints Cll and C12 are the boundary eonditions for the durations of DE and UE 
transmission. 


IV. Resource Allocation Algorithm Design 

The eonsidered optimization problem has a non-eonvex objeetive funetion in fraetional form. 
Besides, the preeoding matriees {Pol^Pul^} and are eoupled in (fT9l) and (l20l) 

leading to a non-eonvex feasible solution set in (l22l) . Also, eonstraints C9 and CIO are eombi- 
natorial eonstraints whieh results in a diseontinuity in the solution set. In general, there is no 
systematie approaeh for solving non-eonvex optimization problems optimally. In many oases, an 
exhaustive searoh method may be needed to obtain the global optimal solution. Yet, applying 
suoh method to our problem will lead to prohibitively high oomputational oomplexity sinoe 
the searoh spaee for the optimal solution grows exponentially with respeot to K and np. In 
order to make the problem traotable, we first transform the objeetive funetion in fraetional form 
into an equivalent objeetive funetion in subtraotive form via fraetional programming theory. 
Subsequently, majorization theory is exploited to obtain the strueture of the optimal preeoding 
poliey to further simplify the problem. Then, we employ eonstraint relaxation to handle the 
binary eonstraints C9 and CIO to obtain an asymptotieally optimal resouree alloeation algorithm 
in high SNR regime and for large numbers of subearriers. 
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A. Transformation of the Optimization Problem 

For notational simplicity, we define IF as the set of feasible solutions of the optimization 
problem in (|2^ spanned by constraints Cl - C12. Without loss of generality, we assume that 
{V,S} G F and the solution set F is non-empty and compact. Then, the maximum energy 
efficiency of the SUDAS assisted communication, denoted as 7]*^, is given by 

(23) 


U{V*,S*) . . U{V,S) 

^eff = = maximize 


Wtp(T’*,iS*) {'p,s}&t Wxp(T’,iS) 

Now, we introduce the following theorem for handling the optimization problem in (l22l) . 

Theorem 1: By nonlinear fractional programming theory [|27ll . the resource allocation policy 
achieves the maximum energy efficiency 17 *^ if and only if it satisfies 


maximize U(V,S) - r]*sUTp{V,S) = U{V\S*) - pl^UppiVfS*) = 0, 

{V,S}&T 

for U{V,S) > 0 and 1{tp{P,S) > 0. 

Proof: Please refer to [l27ll for a proof of Theorem 1. ■ 

Theorem [T] states the necessary and sufficient condition for a resource allocation policy to be 
globally optimal. Hence, for an optimization problem with an objective function in fractional 
form, there exists an equivalent optimization problem with an objective function in subtractive 
form, e.g. U{V*,S*) — r]*gU'Tp{V*,S*) in this paper, such that the same optimal resource 
allocation policy solves both problems. Therefore, without loss of generality, we can focus 
on the objective function in equivalent subtractive form to design a resource allocation policy 
which satisfies Theorem [T] in the sequel. 


B. Asymptotically Optimal Solution 


In this section, we propose an asymptotically optimal iterative algorithm based on the Dinkel- 
bach method ETll for solving (l22l) with an equivalent objective function such that the obtained 
solution satisfies the conditions stated in Theorem [B The proposed iterative algorithm is sum¬ 
marized in Table U (on the next page) and the convergence to the optimal energy efficiency is 
guaranteed if the inner problem (l2^ is solved in each iteration. Please refer to lIZTll for a proof 
of the convergence of the iterative algorithm. 

The iterative algorithm is implemented with a repeated loop. In each iteration in the main 
loop, i.e., lines 3 — 10, we solve the following optimization problem for a given parameter rjes' 


maximize ll{V,S) — rje^UppiV, S) 


s.t. Cl - C12. 


(24) 




TABLE I 


15 


Iterative Resource Allocation Algorithm. 

Algorithm 1 Iterative Resouree Alloeation Algorithm 

1: Initialization: Lmax = the maximum number of iterations and A —0 is the maximum tolerance 
2: Set ?7eff = 0 and iteration index 1 = 0 
3: repeat {Iteration Process: Main Loop} 

4: For a given Tjeff, solve and obtain an intermediate resource allocation policy {V',S'} 

5: if \U{V',S') - rjsaHTp{V',S') < A| then 

6: Convergence =true, return {^*,5*} = {P',5'} and ri*g = g/) 

7: else 

8: Set r/eff = and t = t + 1, convergence = false 

9: end if 

10: until Convergence = true or t = Lmax 


Solution of the Main Loop Problem ll24l) .' The transformed objective function is in subtractive 
form and is parameterized by variable ? 7 eff. Yet, the transformed problem is still a non-convex 
optimization problem. We handle the coupled precoding matrices by studying the structure of the 
optimal precoding matrices for (l24l) . In this context, we define the following matrices to facilitate 
the subsequent presentation. Using singular value decomposition (SVD), the DL two-hop channel 
matrices Hg^g and can be written as 


H 


B-fS 


— iV' 


) ana Lig^gg — u 


A [Lfc] 

S^-UE'^'-S-i-UE 


(V; 


S-)-UE7 


respectively, where Ug_j.g G C 


MxM 


,v|;‘^g e e and v. 


i,k] 

S^UE 


(25) 

g (QMxM 


are unitary matrices. Ag^g and Ag}^gg and are M x N and MxM matrices with main diagonal 


element vectors 


and 


i,k] / [i,k\ 

7s^ue,i V 7s-^UE,2 


M] 
7s^.UE,iJ2 


and 


7b-^s,i Y 7b-^s,2 • • • Y 7b^s,Ri 

respectively, and all other elements equal to zero. Subscript indices Ri = Rank(HgYs 

= Rank(Hg}^gg) denote the rank of matrices Hg_^g and Hg}^gg, respectively. Variables 

\i] \i k] 

7B-ASn 7sYuEn represent the equivalent channel-to-noise ratio (CNR) on spatial channel 
n in subcarrier i of the BS-to-SUDAS channel and the SUDAS-to-UE k channel, respectively. 
Similarly, we can exploit channel reciprocity and apply SVD to the UL two-hop channel matrices 
which yields 


H 


S^B 


= V 


B^-S 


(A 


H 

B^-S 


1 77 


(U; 


*] 

B-:>S 


177 


and H 


i^k] 

UE^S 


= V, 


[i,k] 

S-:>UE 


(A 


i,k\ N 

s->-ue{ 




77 


(26) 


We are now ready to introduce the following theorem. 

Theorem 2: Assuming that Rank(PgL') = Rank(Pgg^) = Rank(Fgg^) = Rank(FgL^) = 
TVs < niin{Rank(Hg}^ug), Rank(Hg^g)}, the optimal linear precoding matrices used at the 
BS and the SUDACs for the maximization problem in (l24l) jointly diagonalize the DL and 
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UL channels of t 
objective function 


le BS-SUDAS-UE link on each subcarrier, despite the non-convexity of the 
1 The optimal precoding matrices have the following structure: 




A 


2,/c] 

Bdl’ 


DL 


_ A I*’''! I'T T 

“ ''S^UE'^'-FDLV^B-s-Sy' 




(27) 


-p[i,A;] 

^UL 


ttW a 

^S-i-UE^'-UEuL’ 


and F 




UL 


- TT^ ’ ‘ A 


(V: 


S^UEy' 1 


(28) 


respectively, where V 3 _^g, and Ug^g are the rightmost columns of Vg^g, Vg)^gg, 

and Ug^g, respectively. Matrices Ag^^^ G G Ag^[ G C^®^^®, and 

A Ml g (J-iVgXiVs 

are diagonal matrices which can be expressed as 




A 


[i,k] 

UEul 


A 


[i,k] 

Ful 


diag (V4^'1 ;e,i • • • 

•^iag (y\lP\]'- e}-,s,i ■ ■ ■ '(/^UE^s.n • • • \lP\]'- e}^s,n^ ; and 

diag(y4% 


(29) 

(30) 

(31) 

(32) 


respectively, where diag(xi,--- ^xk) denotes a diagonal matrix with the diagonal elements 
{xu-- - , xk}. Scalar optimization variables 4^UE.n’ and are, re¬ 

spectively, the equivalent transmit powers of the BS-to-SUDAS link, the SUDAS-to-UE link, 
the UE-to-SUDAS link, and the SUDAS-to-BS link for UE k on spatial channel n and subcarrier 




Proof: Please refer to the Appendix. ■ 

By adopting the optimal precoding matrices provided in Theorem [2l the DE and UE end-to- 
end channel on subcarrier i is converted into parallel spatial channels. More importantly, 
the structure of the optimal precoding matrices simplifies the resource allocation algorithm 
design significantly as the matrix optimization variables can be replaced by equivalent scalar 
optimization variables. As a result, the achievable rates in DE and UE on subcarrier i from the 


H 


^We note that the diagonal structure is also optimal for frequency division duplex systems where 7 ^ and 

gl^g 7 ^ (Hb(^s)^- Only the optimal precoding matrices in l l27b and ( 128b will change accordingly to jointly diagonalize the 


end-to-end channel. 
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BS to UE k via the SUDAS in (fTSl) can be simplified as 


Ns 


DL =£log2 1 
72= 1 

(i + sinrKI) 

, sinrKI 

Ns 

UL 1 

(i + sinrM) 

. sinr“ 


pliM p[i,/ej [i,fcj 

/B-i>S,n-' B^S,n-' S^>UE,n /S^>UE,n 

1 _i_ p[*i^] _i_ p[*i^] ^,[*>^1 

/B^S,n-' B-i.S,n “r -'^S-i^UE.n 7S-l>UE,n 
_[*] p[iM pi^M JiM 

IS^B,n^S^B,n^\JE^S,n /UE^S,n 


i,k] 


.[iM 


1 _|_ _,W p[*:^] I pl'I-^^l 

^ ' IS-^B,n■^S^B,n T“ -'^UE->-S,n7uE^.S,n 


i,k] 


,[iM 


n=l 


(33) 


(34) 


where SINR^lI^ and SINR^’l^ are the received signal-to-interference-plus-noise-ratios (SINRs) 
at UE k and the BS in subcarrier i in spatial subchannel ?7.G{l,...,A7g}, respectively. Although 
the objective function is now a scalar function with respect to the optimization variables, it is still 
non-convex. To obtain a tractable resource allocation algorithm design, we propose the following 
objective function approximation. In particular, the end-to-end DE and UE SINRs on subcarrier 
i for UE k can be approximated, respectively, as 


SINR 


[i,k] 

DLtj, 


SINR 


i^k\ 

DLtt, 


and SINR|5’£^ 


SINR^ 


i,k] 

ULn’ 


sinrS 


\i,k] 


[iM 


7B^-S,n-'B-)-S,n'^S->-UE,n7s-fUE,r^ 


p[^M _i_ p[*i^] 

7B-j-S,n-'^B-i>S,n “T -'^S-)-UE,n7s-fUE,r^ 


SINR 


where 

i,k] 


(35) 


-ii p[*>fc] p[i,^ 

/S-J-B,n-' S^B.n-' UE-i>S,n /UE-i>S,i 


p[*4] [iM 

7S-i>B,n-7S-^>B,r^ “r UE-i>S,n7uE-)-S,; 


We note that this approximation is asymptotically tight for high SNR^ 

The next step is to tackle the non-convexity due to combinatorial constraints C9 and CIO in 

(l22l) . To this end, we adopt the time-sharing relaxation approach. In particular, we relax and 

JiM ; 

'’UL 


in constraints C9 and CIO such that they are non-negative real valued optimization variables 

, i.e., 0 < < a and 0 < s^’l' < /3. It has 


bounded from above by a and f3, respectively 

been shown in ll2^ that the time-sharing relaxation is asymptotically optimal for a sufficiently 

large number of subcarrierj^. Next, we define a set with four auxiliary optimization variables 

'P _ ^p[iM p{iM p{iM p{iM 


B-j-S,n) UE->-S,n> ■' S-)-UE,n 


) Ps^B n} rewrite the transformed objective function in 


as: 


K Tip Nq 

ut,wp,s) = J2'E^ 

k=l i=l n=l 


logs (1 + 


SINR 


■[iM 

'DL„ 


.[iM 


^DL 


+ Sul logs 



(36) 


K np Ns 


-VeS ( PcB+34'EcguDAC+-^-PcuE+y^ ^B^iyS,n+^S^S*4'uE,ri+^fc^UEyS,n+eS^f 

k = l 2=1 T2=l 


i,k] 

S-^B,72 


^It is expected that the high SNR assumption holds for the considered system due the short distance communication between 
the SUDAS and the UEs, i.e., PiEtE,n7KljE.n, » 1- 

^The duality gap due to the time-sharing relaxation is virtually zero for practical numbers of subcarriers, e.g. np > 8 1291 . 
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-[i,k] 


2,/c] 


where SINR^l^ = SINR^l 


UE^S,n 


,\i,k] \i,k] T-i\Uk 

VeUsA , PLkn 


II r T-)\hk] T>H^k] \i,k] r->\i,k] nHM H^k] 

and $ = = ^B^S,n%L > ^s4uE,n = ^s4uE,n%L^ 

H k] \i k] 1 

= ^s4b n'^UL f • auxiliary optimization 


variables in V represent the aetual transmit energy under the time-sharing eondition. As a result, 
the eombinatorial-eonstraint relaxed problem ean be written as: 


s.t. 


Cl: 


C3: 


maximize W-rrans (P,<S) 
v,s 

K np A^s 

k=l i=l n=l 
rip A''s 


K np Nq 


cz^EEE^^ 


< MP 

S^UE,n — -Tmax) 


EE^. 

2=1 n=l 


UE^S,n — -Pmaxfe)V/c, 


k = l 2 = 1 22=1 
K rip Nq 


■^5^] ^ pUL 

S—>^B,r2 — -^max5 


C5-C8, C9: 0 < 4’^’< 


«:EEE^^ 

k = l 2 = 1 22=1 

CIO: 0 < 4l < Cll, C12. (37) 


Optimization problem (1371) is jointly eoneave with respeet to the auxiliary optimization variables 


V and S. We note that by solving optimization problem (iJTl) for -Pg^S’ 


,[i,k] 


,) A 


H,k] 


S^UE,n’ UE-)-S,n’ 

’ and s|;’4 we ean reeover the solution for P^^kn’ 4*-kuE,n’ 4?k-ks,n’ and 
Thus, the solution of (iJTl) is asymptotieally optimal with respeet to (l22l) for high SNR and a 
suffieiently large number of subearriers. 

Now, we propose an algorithm for solving the transformed problem in (1371) . Although the 
transformed problem is jointly eoneave with respeet to the optimization variables and ean be 
solved by eonvex programming solvers, it is diffieult to obtain system design insight from a 
numerieal solution. This motivates us to design an iterative resouree alloeation algorithm whieh 
reveals the strueture of energy-effieient resouree alloeation solutions and serves as a building 
bloek for the suboptimal algorithm proposed in the next seetion. The proposed iterative resouree 
alloeation algorithm is based on alternating optimization. The algorithm is summarized in Table 
Inland is implemented by a repeated loop. In line 2, we first set the iteration index I to zero and 
initialize the resouree alloeation poliey. Variables Pkkkn(0> ^k-iuE,n(0,4*ikn(0,4jEks,n(0) 
(0, ■^urkO’ /^(O denote the resouree alloeation poliey in the l-th iteration. Then, 

in eaeh iteration, we solve (iJTl) . whieh leads to (l38])- (l4^ : 




■ [i,k] jj[i,k] (^[i,k] 

/S-j-UE.n-'^S^-UEp ' “ 


B^-S,) 


H,k] jj[i,k] 
IS^VE,n^S^EE,r, 


Jhk] pM , N 
‘^WB-i>S,n /S->-UE,n-' S^UE,n ^ /B^S,nl 


i,k] 


Q 


[i,k] _ 

B^-S,n“ 


4(1 +W^S[)7BLs,n(l +7KuE,n^k-iuE,J + (7s-iuE,n)^ A + ^efr£B)(Ps^ljE,J^ ln(2) 


\A + \/ln( 2 ) 


(38) 


(39) 
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TABLE II 

Iterative Resource Allocation Algorithm for SUDAS Assisted Communication 


Algorithm 2 Alternating Optimization 


0 


2: Set iteration index i = 0 and initialize (01, ^UB-^s.n(0> and 

/ = i + l 
3: repeat {Loop} 

4: 




Ji.fe] / 








For given -Pg^kuE “ 1) and a{l), solve ( 137b for Pj 
variables 


B-is n using ( 138b which leads to intermediate power allocation 


5: For given „ and a{l), solve for Pt 


k-k^UB rt ''ia equation ( 140b which leads to intermediate power allocation 


variables Pc 




S->UE,n 


6 : Update the DL subcarrier allocation policy via ( 142b with ^{l — 1), PbIJs n’ and a{l) to obtain the intermediate 


,li,k]' 


7: 


8 : 


DL subcarrier allocation policy (Z) 

For given Pg^lfg n “ 1) and /3{l), solve for P, 


i,k] 

UE^S,n 


via equation ( 143 b which leads to intermediate power allocation 


variables Pkg'ls^^ 


9: 


For given Pugi^g „ and P{1), solve ( 137b for Pg_!^B n using ( 144b which leads to intermediate power allocation variables 

p[i,fe]' 

Update the UL subcarrier allocation policy via m with - 1 ), P^'klg.^, and /3{l) to obtain the intermediate 

UL subcarrier allocation policy (Z) 

10: Update a and /3 via standard linear programming methods to obtain intermediate solutions of a' and /3' 

11 : if - 1)1 < N , |P'%L'n - - 1)1 < N , ” 4’l ’ ” 1)1 < 

i4;£s,. - -1)1 < -, i4%L - -1)1 <«, - i)i < ^ 

\a' — a{l — 1)1 < ft, and \P' — p{l — 1)| < ft then 


12 : 

13: 

14: 


Convergence = true, return {Pc 


[i,k]' 


,[i,kY Ji.k]' p[i,kY 


S->UB,rU -f B-s-S.ni *DL > ^ UB->S,.> ^ S->B,n! ^UL 


i,kY 


else 


Convergence 


false. Pc 


[i,k] 


j[i,k]' 


j[i,k] 


S->UB,.(0 — -fs->UB,rtiB-»S,. 

_ plLfc]' „hA]/7\ _ .[i.fc]' 


(Z) = P, 


[i.fc]' 


[i,k] 


B->-S,.> “DL 

= p',i = i + i 

15: end if 

16. until Z = Pmax 

,[Lfc] 


(Z) = s 


[i,kY p[i,fc] 

DL 5 -* UE->S,n 


(0 = 


with .Ps%UEn(^~l) from the last iteration, where [x]’*' = max{x, 0}. Then, the obtained 


\i k] 

intermediate power alloeation variable „ is used as an input for solving 
via the following equations: 

2 


for Pc 


[i,k] 


S^UE.n ~ 


p[iM 


o[*,fc] 

“S-tUE,i 


Ji] 

/B-tS.n-'^B^S,: 


I ) 

/B-tS,n /S^UE,n-'^B-tS,n ^ /S^UE.n) 


f.Zc] 


{i,k] 


S^UE.n 


(40) 


Q 


i,k] 

S^UE,n ' 


(TBkj,S,n )^(Z+ i?el£s)(fS}'s,„)Un(2) + ( 7 . 


pl;l„ + l)4(l+ufi)7K’uE.. 


[Zc] \ 


B—tS.n'' B—tS,n 


(41) 


Vh + r}efr£s\/ln(2) 

-dUl) are obtained by standard convex optimization techniques. A and 5 in (IMl) and (l40l) 
are the Lagrange multipliers for constraints Cl and C2 in (lT7l) . respectively. In particular, A and 
6 are monotonically decreasing with respect to Pg^^s n -^s-^ue n’ respectively, and control 












20 


the transmit power at the BS and the SUDAS to satisfy eonstraints Cl and C2, respeetively. 
Besides, > 0 is the Lagrange multiplier assoeiated with the minimum required DL data rate 
eonstraint C5 for delay sensitive UE k. The optimal values of A, 5, and iteration 

ean be found by a standard gradient algorithm sueh that eonstraints Cl, C2, and C4 in (iTTl) are 
satisfied. Variable r/eg > 0 generated by the Dinkelbaeh method prevents unneeessary energy 
expenditures by redueing the values of ^s^uE,n and (l4T]) . respeetively. 

Besides, the power alloeation strategy in (l3^ and (l40l) is analogous to the water-filling solution 
in traditional single-hop eommunieation systems. In partieular, n n water 

levels for eontrolling the alloeated power. Interestingly, the water level in the power alloeation 
for the BS-to-SUDAS link depends on the assoeiated ehannel gain whieh is different from the 
power alloeation in non-SUDAS assisted eommunieation ffTSll . [fT^ . Furthermore, it ean be seen 
from (l39l) and (1411) that the water levels of different users ean be different. Speeifieally, if the 
end-to-end ehannel gains of two users are the same, to satisfy the data rate requirement, the 
water level of the delay-sensitive user is generally higher than that of the non-delay sensitive 
user. 

After obtaining the intermediate DL power alloeation poliey, ef. lines 4, 5, we update the DL 
subearrier alloeation, of. line 6, as: 

a if k = arg(1 + (En=ilog 2 (l + - 

0 otherwise 




sinr^’l: 




.1421 



Here, SINR^l^ is obtained by substituting the intermediate solution of „ and L’g^^uE 
i.e., (l38l) and (l40l) . into (1351) in the l-th iteration. We note that the optimal value of of the 
relaxed problem is a diserete value, of. (I42l) . i.e., the eonstraint relaxation is tight. 

H k] \i k] 

Similarly, we optimize the UL power alloeation variables, n and Tg2^B n’ sequentially, 

of. lines 7, 8, via the following equations: 


p[®4] _ 

-'^UE-i>S,n “ 




PA _ [i] p[i,fc] _ „ 

/S-i>B,n-' S^B,n I “UE^.S,n 7S-i>B,n-'^S-i>B,n ^ 

n( [A [Ak] p[i,k] , [i,k] ^ 

PlS^B,nlVE^S,n-^S^B,n 7uE^-S,nJ J 

■ li,k] p[i,fc] f(^[i,k] _Ai,k] Pj-M _rf 
7uE^S,n-'^UE^S,n l“S^>B,n 7uE-i>S,n-' UE^S,n ^ 


jAk] pg,Kj , m N 
^ WS-)-B,n7uE-)-S,n-^UE-)-S,n ^ 7s^.B,nl 


i,k] 


(43) 

(44) 
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respectively, where 


o[*,fc] 

“UE^.S,n 


(7| 




'(l/lfc + r7eff£fc)(-Pc 


i,k] 


ln(2) + ( 7 . 


-p[i,k] 

S^-B,n-'^S-i>B,; 


+1)4(1 + ^r'uL)7uE^s, 


^(45) 


\/7fe + ^eff£fc\/ln(2) 


“S-)-B,; 


4(1 + ^l;D7^^B,n(l + 7|;E^S,n^/;E^S. J + (7 |;e^S, 1^(2) 

V'^ + 77effes\/ln(2) 


(46) 


'ijjk and 0 in (l45l) and (l46l) are the Lagrange multipliers with respect to power consumption 
constraints C3 and C4 in (l37l) . respectively. Besides, is the Lagrange multiplier associated 
with the minimum required UL data rate constraint C6 for delay sensitive UE k. The optimal 
values of 0^, 0, and te|jL in each iteration can be easily obtained again with a standard gradient 
algorithm such that constraints C3, C4, and C6 in (l37l) are satisfied. 

\i k] 

Then, we update the UL subcarrier allocation policy Sul via 

lik] \ (3 if k = arg max (1 + w^’l) (En=ilog 2 (l + SINRul, 

I 0 otherwise 


SINK 


i+sinrSK 


.1471 


where SINRul^ is obtained by substituting the intermediate solutions for n -^u&Ls n’ 
i.e., (l43l) and (l44l) . into (l35l) in the l-th iteration. Again, the constraint relaxation is tight. 

Subsequently, for a given UL and DL power allocation policy and given and the 
optimization problem is a linear programming with respect to a and 0. Thus, we can update 
a and (5 via standard linear programming methods to obtain intermediate solutions for a' and 
13'. Then, the overall procedure is repeated iteratively until we reach the maximum number of 
iterations or convergence is achieved. We note that for a sufficient number of iterations, the 
convergence to the optimal solution of (l37l) is guaranteed since (l37l) is jointly concave with 
respect to the optimization variables ll30]l . Besides, the proposed algorithm has a polynomial 
time computational complexity. 


C. Suboptimal Solution 

In the last section, we proposed an asymptotically globally optimal algorithm based on the high 
SNR assumption. In this section, we propose a suboptimal resource allocation algorithm which 
achieves a local optimal solution of (l22l) for arbitrary SNR values. Similar to the asymptotically 
optimal solution, we apply Theorem [U Theorem |2l and Algorithm [T] to simplify the power 
allocation and subcarrier allocation. In particular, the DL and UL data rates of UE k on subcarrier 
i are given by (l3^ and (l34l) . respectively. It can be observed that (l3^ and (l34l) are concave 
functions with respect to -^s*^UEn’ -^u&Lsn individually, when the other 
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variables are fixed. Thus, we ean apply alternating optimization to obtain a loeal optimal solution 
ifMI of (l22l) . We note that unlike the proposed asymptotieally optimal seheme, the high SNR 
assumption is not required to eonvexify the problem. The suboptimal solution ean be obtained by 
Algorithm [2] in Table UH but now, we update the power alloeation variables, i.e., lines 4,5,7,8, 
in Algorithm 2, by using the following equations: 
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whieh are obtained by applying standard optimization teehnique lf24]| . Besides, the subearrier 
alloeation polieies for DL and UL are still given by (l42l) and (l47l) . respeetively, exeept that 

, respeetively. The optimization 
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we replaee SINRI^l; , SINR;;l; by SINR^J^, SINR^^’l^ in 
variables are updated repeatedly until eonvergenee or the maximum number of iterations is 
reaehed. In eontrast to the asymptotieally optimal algorithm in Seetion IIV-BI whieh may not 
even aehieve a loeally optimal solution for finite SNRs, the suboptimal iterative algorithm is 
guaranteed to eonverge to a loeal optimum If^ for arbitrary SNR values. 

V. Results and Discussion 

In this seetion, we evaluate the system performanee based on Monte Carlo simulations. We 
assume an indoor environment with AT = 4 UEs and M SUDACs and an outdoor BS. The 
distanees between the BS and UEs and between eaeh SUDAS and eaeh UE are 100 meters 
and 4 meters, respeetively. Eor the BS-to-SUDAS links, we adopt the Urban maero outdoor-to- 
indoor seenario of the Wireless World Initiative New Radio (WINNER+) ehannel model ifSTll . The 
eenter frequeney and the bandwidth of the lieensed band are 800 MHz and 20 MHz, respeetively. 
There are np = 1200 subearriers with 15 kHz subearrier bandwidth resulting in 18 MHz signal 
bandwidth for data transmissioiu. Henee, the BS-to-SUDAS link eonfiguration is in aeeordanee 


[i,k] 


^The proposed SUDAS can be easily extended to the case when carrier aggregation is implemented at the BS to create a 
large signal bandwidth 100 MHz) in the licensed band. 
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with the system parameters adopted in the Long Term Evolution (LTE) standard As for 
the SUDAS-to-UE links, we adopt the IEEE 802.Had ehannel model [[^ in the range of 60 
GHz and assume that M orthogonal sub-bands are available. The maximum transmit power of 
the SUDACs and UEs is set to MP^ax = -Pmax = Pmaxk = 23 dBm whieh is in accordance 
with the maximum power spectral density suggested by the Harmonized European Standard for 
the mmW frequency band, i.e., 13 dBm-per-MHz, and the typical maximum transmit power 
budgets of UEs. Eor simplicity, we assume that Ns = min{A^, M} for studying the system 
performance. We model the SUDAS-to-BS and UE-to-SUDAS links as the conjugate transpose of 
the BS-to-SUDAS and SUDAS-to-UE links, respectively. Also, the power amplifier efficiencies 
of all power amplifiers are set to 25%. The circuit power consumption for the BS and each 
BS antenna are given by Fcb = 15 Watt and PAnts = 0.975 Watt, respectively ll34]l . [|35l. 
The circuit power consumption per SUDAC and UE are set to Pcsudac = 0.1 Watt ll^ and 
PcuB = 1 Watt, respectively. We assume that there is always one delay sensitive UE requiring 
=20 Mbit/s and =20 Mbit/s in DE and UE, respectively. Also, Ns is chosen 

as Ns = min{Rank(Hg^ug),Rank(HB^g)}. All results were averaged over 10000 different 
multipath fading channel realizations. 

A. Convergence of the Proposed Iterative Algorithm 

Eigure |3] illustrates the convergence of the proposed optimal and suboptimal algorithms for 
N = 8 antennas at the BS and M = 8 SUDACs, and different maximum transmit powers at the 
BS, Pt- We compare the system performance of the proposed algorithms with a performance 
upper bound which is obtained by computing the optimal objective value in (iTTl) for noise-free 
reception at the UEs and the BS. The performance gap between the asymptotically optimal 
performance and the upper bound constitutes an upper bound on the performance loss due 
to the high SINK approximation adopted in (1351) . The number of iterations is defined as the 
aggregate number of iterations required by Algorithms 1 and 2. It can be observed that the 
proposed asymptotically optimal algorithm approaches 99% of the upper bound value after 
20 iterations which confirms the practicality of the proposed iterative algorithm. Besides, the 
suboptimal resource allocation algorithm, achieves 90% of the upper bound value in the low 
transmit power regime, i.e., Rp = 19 dBm, and virtually the same energy efficiency as the upper 
bound performance in the high transmit power regime, i.e., Fp = 46 dBm. In the following 
case studies, the number of iterations is set to 30 in order to illustrate the performance of the 
proposed algorithms. 
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Fig. 3. Average energy efficiency (Mbits/Joule) versus 
the number of iterations for different maximum transmit 
power budgets at the BS. 



Fig. 4. Average energy efficiency (Mbits/Joule) versus 
the maximum transmit power at the BS (dBm) for different 
communication systems. The double-sided arrows indicate 
the performance gain achieved by the proposed SUDAS. 


B. Average System Energy Efficiency versus Maximum Transmit Power 


Figure |4] illustrates the average system energy efficieney versus the maximum DL transmit 
power at the BS for M = 8 SUDACs for different systems and = 8 BS antennas. It ean 
be observed that the average system energy effieieney of the two proposed resource allocation 
algorithms for SUDAS is a monotonically non-decreasing function of Ft- In particular, starting 
from a small value of Fp, the energy efficiency increases slowly with increasing Fp and 
then saturates when Fp > 37 dBm. This is due to the fact that the two proposed algorithms 
strike a balance between system energy efficiency and power consumption. In fact, once the 
maximum energy efficiency of the SUDAS is achieved, even if there is more power available 
for transmission, the BS will not consume extra DL transmit power for improving the data 
rate, cf. (1381) . This is because a further increase in the BS transmit power would only result in a 
degradation of the energy efficiency. Moreover, we compare the energy efficiency of the proposed 
SUDAS with a benchmark MIMO system and a baseline system. We focus on two system design 
objectives for the reference systems, namely, system throughput maximization (TP-Max) and 
energy efficiency maximization (EE-Max). For the benchmark MIMO system, we assume that 
each UE is equipped with N receive antennas but the SUDAS is not used and optimal resource 
allocation is performecj^ Besides, we assume that the circuit power consumption at the UE does 
not scale with the number of antennas and only the licensed band is used for the benchmark 


®The optimal resource allocation for the benchmark system can be obtained by following a similar method as the one proposed 
in this paper applying also fractional programming and majorization theory. 
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system. In other words, the average system energy efficieney of the benchmark system serves 
as a performance upper bound for the proposed SUDAS. For the baseline system, we assume 
that the BS and the single-antenna UEs perform optimal resource allocation and utilize only the 
licensed frequency band, i.e., the SUDAS is not used. As can be observed from Figure H for 
high BS transmit power budgets, the SUDAS achieves more than 80% of the benchmark MIMO 
system performance even though the UEs are only equipped with single antennas. Also, the 
SUDAS provides a huge system performance gain compared to the baseline system which does 
not employ SUDAS since the proposed SUDAS allows the single-antenna UEs to exploit spatial 
and frequency multiplexing gains. On the other hand, increasing the number of BS antennas 
dramatically in the baseline system from = 8 to = 64, i.e., to a large-scale antenna system, 
does not necessarily improve the system energy efficiency. In fact, in the baseline system, the 
higher power consumption, which increases linearly with the number of BS antennas, outweighs 
the system throughput gain, which only scales logarithmically with the additional BS antennas. 

Eigure [5] depicts the average time allocation for DE and UE transmission. It can be observed 
that the optimal time allocation depends on the transmit power budget of the systems. In 
particular, when the power budget of the BS for DE communication is small compared to 
the total transmit power budget for UE communication, e.g. Pt < 28 dBm, the period of time 
allocated for DE transmission is shorter than that allocated for UE transmission. Because of the 
limited power budget and the circuit power consumption, it is preferable for the BS to transmit 
a sufficiently large power over a short period of time rather than a small power over a longer 
time to maximize the system energy efficiency and to fulfill the data rate requirement of the DE 
delay sensitive UEs. On the contrary, when the power budget of the BS is large compared to 
that of the UEs, the system allocates more time resources to the DE compared to the UE, since 
the BS can now transmit a large enough power to compensate the circuit power consumption 
for a longer time span to maximize the system energy efficiency. 

C. Average System Throughput versus Maximum Transmit Power 

Eigure [6] illustrates the average system throughput versus the maximum transmit power at 
the BS for A^ = 8 BS antennas, Ff = 4 UEs, and M = 8. We compare the two proposed 
algorithms with the two aforementioned reference systems. The proposed SUDAS performs 
closely to the benchmark scheme in the low DE transmit power budget regime, e.g. Pt < 31 dBm. 
This is due to the fact that the proposed SUDAS allows the single-antenna UEs to transmit or 
receive multiple parallel data streams by utilizing the large bandwidth available in the unlicensed 
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Fig. 5. Average DL and UL transmission durations versus 
the maximum transmit power at the BS (dBm). 


Fig. 6. Average system throughput (Mbits/s) versus the 
maximum transmit power at the BS (dBm). 


band. Besides, for all eonsidered systems, the average system throughput of all the systems 
inereases monotonieally with the maximum DL transmit power Pt- Yet, for the systems aiming 
at maximizing energy effieieney, the eorresponding system throughput saturates in the high 
transmit power allowanee regime, i.e., Pp > 37 dBm, in the eonsidered system setting. In faet, 
the energy-effieient SUDAS does not further inerease the DL transmit power sinee the system 
throughput gain due to a higher transmit power eannot eompensate for the inereased transmit 
power, i.e., the energy effieieney would deerease. As for the benehmark and baseline systems 
aiming at system throughput maximization, the average system throughput inereases with the DL 
transmit power without saturation. For system throughput maximization, the BS always utilizes 
the entire available DL power budget. Yet, the inereased system throughput eomes at the expense 
of a severely degraded system energy effieieney, ef. Figure IH 

D. Average System Performance versus Number of SUDACs 

Figures |7] and [8] illustrate the average energy effieieney and throughput versus the number of 
SUDACs for = 8 BS antennas, Pp = 37 dBm, and different systems. It ean be observed 
that both the system energy effieieney and the system throughput of the proposed SUDAS 
grows with the number of SUDACs, despite the inereased power eonsumption assoeiated with 
eaeh additional SUDAC. For N > M, for DL transmission, additional SUDACs faeilitate a 
more effieient eonversion of the spatial multiplexing gain in the lieensed band to a frequeney 
multiplexing gain in the unlieensed band whieh leads to a signifieant data rate improvement. 
Similarly, for UL transmission, the SUDACs help in eonverting the frequeney multiplexing gain 
in the unlieensed band to a spatial multiplexing gain in the lieensed band. For M > N, inereasing 
the number of SUDACs in the system leads to more spatial diversity whieh also improves energy 
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Fig. 7. Average energy efficiency (Mbits/Joule) versus Fig. 8. Average system throughput (Mbits/s) versus the 

the number of SUDACs. number of SUDACs. 


efficiency and system throughput. Besides, a substantial performance gain can be achieved by 
the SUDAS compared to the baseline system for an increasing number of available SUDACs. 

VI. Conclusions 

In this paper, we studied the resource allocation algorithm design for SUDAS assisted outdoor- 
to-indoor communication. Specifically, the proposed SUDAS simultaneously utili z es licensed 
and unlicensed frequency bands to facilitate spatial and frequency multiplexing gains for single¬ 
antenna UEs in DL and UL, respectively. The resource allocation algorithm design was formu¬ 
lated as a non-convex matrix optimization problem. In order to obtain a tractable solution, we 
revealed the structure of the optimal precoding matrices such that the problem could be trans¬ 
formed into a scalar optimization problem. Based on this result, we proposed an asymptotically 
globally optimal and a suboptimal iterative resource allocation algorithm to solve the problem 
by alternating optimization. Our simulation results showed that the proposed SUDAS assisted 
transmission provides substantial energy efficiency and throughput gains compared to baseline 
systems which utilize only the licensed frequency spectrum for communication. 

Appendix-Proof of Theorem [2] 

Due to the page limitation, we provide only a sketch of the proof which follows a similar 
approach as in [l25l, jlTTll and uses majorization theory. We show that the optimal precoding and 
post-processing matrices jointly diagonalize the DL and UL end-to-end channel matrices on each 
subcarrier for the maximization of the transformed objective function in subtractive form in (l24l) . 
Lirst, we consider the objective function in subtractive form for UL k on a per-subcarrier basis 
with respect to the optimization variables. In particular, the per-subcarrier objective function for 
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UE k consists of two parts, 


hiv.s) 

f2{V,S) 


= -sfen logs ( det[E|^’L^]) - logs ( det[E|^’L^]) 

= Tr + 4 l Tr (g&l") + 4l 

+ £,4"lTr(p&^l(P&^l)^), 


egTr 



( 52 ) 


(53) 


such that the maximization of the per subearrier objeetive funetion ean be expressed as 

minimize - fi {V, S) + rjes f 2 {,'P, S). (54) 


Besides, the determinant of the MSE matrix on subearrier i for UE k ean be written as 

Ns 
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DL 
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where [X]a^6 extraets the (a,6)-th element of matrix X. fi{V,S) is a Sehur-eoneave funetion 
with respeet to the optimal preeoding matriees [l25]l for a given subcarrier alloeation poliey S. 
Thus, —fi{V,S) is minimized when the MSE matrix E^’lI is a diagonal matrix. Eurthermore, 
the traee operator in f 2 {V,S) for the eomputation of the total power eonsumption is also 
a Sehur-eoneave funetion with respeet to the optimal preeoding matriees. Thus, the optimal 
preeoding matriees for the minimization of funetion f 2 {V,S) should diagonalize the input 
matrix of the traee funetion, of. [|^ Chapter 9.B.1] and [|^ Chapter 9.H.l.h]. Similarly, the 
power eonsumption funotions on the left hand side of oonstraints C1-C4 in (l24l) are also Sehur- 
eoneave funotions and are minimized if the input matriees of the traee funotions are diagonal. 
Besides, the non-negative weighted sum of Sehur-eoneave funotions over the subearrier and 
UE indices preserves Sohur-oonoavity. In other words, the optimal preeoding matriees should 
jointly diagonalize the subtraetive form objeetive funetion in (l24l) and simultaneously diagonalize 
matriees Pul^(Pul^)^, Pbl'(Pdl')^) and This observation establishes a neeessary 

oondition for the struoture of the optimal preeoding matriees. Einally, by performing SVD on the 
ehannel matriees and after some mathematieal manipulations, it ean be verified that the matriees 
in (1271) and (1281) satisfy the optimality oondition. 
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